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ABSTRACT 

Be stars possess gaseous circumstellar disks that modify in many ways the spectrum of the central B star. 
Furthermore, they exhibit variability at several timescales and for a large number of observables. Putting the 
pieces together of this dynamical behavior is not an easy task and requires a detailed understanding of the phys- 
ical processes that control the temporal evolution of the observables. There is an increasing body of evidence 
that suggests that Be disks are well described by standard ff-disk theory. This paper is the first of a series that 
aims at studying the possibility of inferring several disk and stellar parameters through the follow-up of vari- 
ous observables. Here we study the temporal evolution of the disk density for different dynamical scenarios, 
including the disk build-up as a result of a long and steady mass injection from the star, the disk dissipation that 
occurs after mass injection is turned off, as well as scenarios in which active periods are followed by periods of 
quiescence. For those scenarios, we investigate the temporal evolution of continuum photometric observables 
using a 3-D non-LTE radiative transfer code. We show that lightcurves for different wavelengths are specific of 
a mass loss history, inclination angle and a viscosity parameter. The diagnostic potential of those lightcurves 
is also discussed. 

Subject headings: circumstellar matter radiative transfer stars: emission-line, Be 



1. INTRODUCTION 

Line emission in most stellar spectra arises from ionized 
gas beyond the photospheric level. In particular, emission-line 
spectra of B-type stars (e.g., Be, B[e], Herbig Ae/Be) are due 
to extended circumstellar (CS) envelopes. Additionally, the 
ionized CS gas produces co ntinuum radiation due to free-fre e 
and free-bound transitions dGehrz. Hackwell. & Jonesll 19741) . 
The envelope contribution to the object's brightness depends 
on the distribution of the density, temperature, and ionization 
degree throughout the CS envelope. Classical Be stars are 
a large class of objects in which CS contribution to the stel- 
lar continuum can be significant. They can display strong IR 
excesses, which can be up to « 40 times larg er than the pho- 
tosph eric flux of the central star at 12 fxm dBeichman et alj 
1988). 

In the past decade or so, a consensus has emerged that the 
CS matter around Be stars is distributed in a disk. The pres- 
ence of rotating, flattened material (i.e. a disk) was initially 
inferred from the typical double-peaked emission line pro- 
files seen in many Be stars and has since been confirmed 
by modern high-angular resolution techniques which have 
resolved the disks arou nd several nearby Be stars (see e.g. 
lOuirrenbach et al.ll 19971) . 

The origin of those disks has been the subject of much de- 
bate. Clearly, a necessary ingredient for disk formation is the 
star's rapid rotation rate, which allows material to be more 
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easily lifted off the surface of the star. Nevertheless, since 
mos t Be stars seem not to be rotating at their break-up speed 
(e.g. lCranmerl2005b, other mechan ism(s) are likely necessary 
(for an overview see lOwoc ki 2006). Regardless of how mate- 
rial is ejected into orbit, another mechanism is required to dis- 
tribute the material throughout the disk. The viscous decr etion 
disk model (VDDM) , first sugge sted bv lLee et alj d 199 lb and 
furt her developed by iBTorkmarJ (11997b . lOkazaki et alj (12002b 
and lBiorkman & Carciofil (12005b . among others, uses the an- 
gular momentum transport by turbulent viscosity to lift ma- 
terial into higher orbits, thereby causing the disk to grow in 
size. A theoretical prediction of this model is that material is 
in Keplerian rotation throughout the disk. Spectrointerferom- 
etry and spectroastrometry are starting to provide clear evi- 
den ce that, in most obser ved and so far analyzed systems (see 
e.g. Meilland et al. 201 1 ), the disks rotate in a Keplerian fash- 
ion dMeilland et alj |2007l: IStefl et alj|20!lt iKraus et alj|20!lt 
iDelaa et alJl201 lb This important fact, together with other 
observational facts outlined in Carciofi 2011, are properties 
that only the VDDM can reproduce. 

Photometric observations o ffer a possibility to s tudy differ- 
ent regions of Be disks (e.g. ICarciofi et alj 120061) . At short 
wavelengths, say, V-band, excess continuum radiation arises 
from a relatively small area near the star, whereas for long 
wavelengths the emission area increase s with the wavelength 
appro ximately as /I 8 ' 11 (see Eq. A15 of ICarciofi & Bjorkmanl 
12006b . Some stars are known to have had a nearly stable con- 
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tinuum emission for decades. C Tauri being a notable example 
dStefl et al.ll2009l) . This suggests that these stars possess a de- 
cretion disk that has been fed at a roughly steady rate. How- 
ever, there are also cases that a star that has been stable for a 
long time suddenly looses its visible and IR excesses (e.g., n 
Aqr, IWisniewsk i et al. 2010). This is interpreted as the dis- 
sipation of the disk by some mechanism, as a result of the 
mass loss from the star being turned off. Similarly, there are 
cases where the disk has been rebuilt after dissipation, and is 
now stable again. During the time of disk growth/dissipation 
there is often short term, small scale photometric variabil- 
ity, a result of transitory outbursts with timescales of days 
or weeks (sometimes c alled "flickering activity", as in Cen, 
lHanuschik et al.ll993l) . Finally, a good fraction of Be stars are 
intrinsically variable at different timescales. Their photomet- 
ric variability can be found periodic, quasi-periodic or irregu- 
lar dMennickent et al.ll994tlsterken et al.ll996HSabogal et al.l 



|2008 )_ or can exhibit episodic outbursts of variable dura- 
tion ([Hubert et al.H2000UMennickent et alj|2002l) . Recently, a 
time-dependent model akin to the one employed in this work 
was successfull y applied for modeli ng the dissipation of the 
disk of 28 CMa (ICarciofi et alj|2012l) . From that analysis, the 
authors were able to determine the viscosity parameter of the 
disk, having found a — 1 .0 ± 0.2. 

The above suggests the existence of at least two timescales 
controlling the disk photometric properties: 1) T m , timescale 
for the variability of the mass injection into the disk, related 
to the rate of stellar mass ejection events and the length of 
these events, and, 2) t^, timescale for the disk to redistribute 
the injected material. 

To critically test the VDDM against observations, the 
structure of the disk must be determined from hydrody- 
namical equations. Using constant mass injection rates3 
and non-LTE codes based on different approa ches (e.g. 
ICarciofi & Bjorkmanl l2006t ISigut & Jonel 120071) . the first 
quantitative tests of the near stea dy state solution ha ve been 



successfully carried out on 5 Sco (ICarciofi et al 
dTvcner et all 120081) . and £ Tau ( ICarciofi et al.1 



12006b . yQph 
20091) . How- 



ever, those models can only be applied to objects that went 
through a sufficiently long and stable decretion phase. 

The purpose of this paper is to provide a framework to un- 
derstand the effects of time variable mass loss rates on the 
structure of the disk and its observational consequences. To 
accomplish this we will first describe the model in § 12 then 
in §[3]we study two limiting cases, for which the timescale of 
mass loss rate variability is either much longer or comparable 
to the disk timescale introduced above. In particular, we study 
the growth of a disk where none has been previously present 
and the dissipation of a pre-existing disk. In the intermediate 
regime where those timescales are comparable we study sin- 
gle outburst-like events and periodic ones. In §|4]the resulting 
photometric observables are described, which are then com- 
pared to previous work and observations in § [5] We finally 
summarize this work in § [6] 

2. MODEL DESCRIPTION 

The tempo r al evo l ution of Be di s ks was studied by 
Okazaki et all (I2002T) . lOkazakil (120071) Uones et al.l d2008l) . 
Sigut & Jonesl (120071) and ICarciofi etal.1 (120121) . With the ex- 



ception of the former, these studies focused on the evolution 
of the surface density of a disk fed by a constant mass injec- 

1 In this paper, we distinguish the mass loss rate and the mass injection 
rate which is the normalized rate at which mass is injected into the disk. 



tion rate and constrained the power-law index n describing the 
radial dependency of the density (p(r) oc r~ n ). In isothermal 
disks, lOkazakil (120071) found that n is always larger than 7/2, 
but approaches this value as time goes to infinity. In order 
to analyze the observational signatures of variable mass in- 
jection ra tes, we use the time-dependent hyd rodynamic code 
singlebe (IOkazakill2007T:IOkazaki et al.ll2002l) . This code per- 
forms one-dimensional (1-D) simulations of the structure and 
evolution of an isothermal viscous de cretion disk by solving 
the ti me-dependent fluid equations dLynden-Bell & Pringld 
119741) in the thin disk approximation. The evolution of such a 
disk is described by the following one-dimensi onal, diffusion - 
type equation of the surface density 2 (see e.g.. lPringlell981l) . 



di 



ld_ 

r dr 



(1) 



where a is the Shakura-Sunyaev viscosity parameter, c s is the 
isothermal sound speed, and Q is the angular frequency of the 
disk rotation. In this equation, we take the angular frequency 
of disk rotation to be circularly Keplerian, i.e., Q(r) = Q.^, 
where D. K = (GM/c 3 ) 1 ' 2 - 

We adopt the (torque-free) outflow boundary condition, 
£ = 0, for both the inner and outer edges of the disk. We 
further assume that mass is injected at a point located just 
above the photosphere, r,„ = 17?*. We place the inner bound- 
ary inside the surface of the star, while the outer boundary is 
placed at r out = 1000/?*, which is roughly the location of the 
trans-sonic point for photo-evaporation of the disk. Note that 
for steady-state outflow, the surface density of the disk scales 
as S oc (Ml a) V>w7#I (see eq. 37 lBiorkmanlll997l) . Conse- 
quently, if a different physical mechanism is responsible for 
truncating the disk (for example, a binary companion), our re- 
sults for the asymptotic values of the surface density should 
be scaled by y o «r/1000fl*. 

The output of singlebe is the surface density, S(r, f), as a 
function of radius and time for a given stellar mass loss history 
and a viscosity parameter. This surface density is converted 
to volume density using the usual vertical hydrostatic equi- 
librium solution (a Gaussian) with a power law scale height, 
H = Hoir/R*) 15 . The volume density is then used as input for 
the three-dime nsional non-LTE Monte Carlo rad iative trans- 
fer code HDUST (ICarciofi & Biorkmanll2006ll2008l) . This code 
performs a full spectral synthesis by simultaneously solving 
the NLTE statistical equilibrium equations and radiative equi- 
librium equation to obtain the hydrogen level populations and 
electron temperature throughout the disk. The code's output 
consists of the SED, polarization spectrum and other observ- 
ables of interest. 

In this work, we adopt a rotationally deformed and gravity 
darkened star whose parameters are typical of a B2Ve star 
(Table Q]). To account for the effects of rotation, the star is 
divided in a number of latitude bins (typically 100), each with 
its effective temperature, gravity and a spec tral shape give n by 
the appropriate Kurucz model atmosphere dKuruczll 19941) . 

To specify the disk structure, we need, as input for singlebe, 
the mass loss rate as a function of time, M(t), and a. Here we 
define the stellar mass loss rate as the mass injected per unit of 
time at the stellar radius. In order to be able to compare differ- 
ent mass loss rate scenarios, we introduce the parameter Zo, 
which is the density at the base of the disk after a sufficient 
long decretion phase. In this case, th e surface density ap- 
proaches the near steady state solution (iBiorkman & Carciofil 
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TABLE 1 

Stellar main parameters used in the 
seviulations. 



Parameter 


Value 


Mass 


9.0 M Q 


Polar radius S po ic 


5.7 R Q 


Equatorial radius 


R< = 6.5R a 


Rotation speed 


113 km/s 


Keplerian speed at equator 


V K = 514 km/s 


Breakup speed ^J^GM/R po i 


V c = 448 km/s 


a /sic 


0.8 


Oblateness 


1.14 


Polar temperature 


22000 K 


Luminosity 


5980 Lq 


TABLE 2 




Disk parameters used in the simulations. 


Parameter 


Value 


So, surface density at the base of the disk 


0.85 g.cirT 2 


M [KT 9 M /year] 


1.64, 4.90, 8.15, 11.41 and 16.29 


a 


0.1,0.3, 0.5, 0.7 and 1.0 


Flaring parameter 


1.5 



120051) . 

E(r) = ^(rlR+T 2 . (2) 

Because Eo depends on the ratio M/a, these quantities are not 
independent parameters in our models. Hence, given a value 
for E , a value of a is uniquely associated with a value of M 
and vice-versa. 

Table [2] summarizes the disk parameters adopted in this 
work. We chose the value Eo = 0.85 g cirT 2 , which corre- 
sponds to a volume density of 3 xlO~ 11 g cm 3 , a typical valu e 
for dense Be disks (ICarciofi et al.ll2006t IWaters et al.||1988l) . 
Several values of a were investigated in the range 0.1 — 1 .0, 
whose corresponding mass loss rates lie in the range 1.6 — 
16 x 10~ 9 M yr~' . We specify here that a is constant through 
all the radii. The possibility of a radial dependence of a is 
mentioned in Sect. 14.21 

3. TEMPORAL EVOLUTION OF THE SURFACE DENSITY 

Our purpose is to guide the analysis of Be disk observa- 
tions by exploring simple dynamical models and their reper- 
cussions on the photometric observables. As discussed in the 
introduction, observations indicate that the state of the disk, 
and hence its photometric properties, is controlled by two 
competing timescales, T; n and tj. This tells us that there will 
be three regimes: 

1. t;„ ^ t&- This situation is irrelevant for the purpose of 
this paper because it does not produce significant ob- 
servable effects. 

2. Ti n » Td. Here, as limiting cases we study the creation 
of a new disk fed at a constant mass injection rate and 
the dissipation of a pre-existing disk. 

3. Tin ~ T&- In this case there will be an interesting inter- 
action between the two competing timescales. In this 
work we explore periodic mass injection scenarios and 
episodic ones. 



3.1. Reference cases: disk growth and dissipation 
3.1.1. Disk growth 

A decretion dis k never experie nces steady state: it either 
grows or decays (|Okazaki 2007). However, it is useful to 
study the solution of a disk subject to a constant mass injec- 
tion rate when time goes to infinity. Even though steady state 
is never physically realized, a limit value for the surface den- 
sity of an ideally grown, unperturbed and isothermal disk can 
be determined analytically, in w hich case the surface densit y 
takes the simple form of Eq. (O dBjorkman & Carciofill2005l) . 

Figure [TJ left panels, shows the temporal evolution of the 
surface density of a newly formed disk subject to a constant 
injection rate for 50 years. The steady-state limit of Eq. (fJJ is 
shown as the thick line. In the entire disk the surface density 
grows with time, albeit with very different rates depending 
on the distance from the star. Therefore, the time required to 
reach quasi-steady state depends strongly on the location in 
the disk as well as on the viscosity parameter a. Also, the 
slope of the surface density profile is very steep at the early 
stages of disk formation and then slowly decreases with time, 
asymptotically reaching the steady-state slope. 

In this simple situation of a disk fed at a constant rate, the 
influence of the viscosity parameter a on the disk is simply 
to scale time up and down. It is immediately evident from 
Eq. ([TJ that changing a is equivalent to changing the time, 
as long as the ratio between the orbital speed and the sound 
speed is kept constant. Therefore, the higher the a the faster 
the disk will grow. For instance, a disk with a — 1 will evolve 
strictly 10 times faster than a disk with a — 0.1. 

Naturally, one could expect that the disk formation occurs 
in the viscous diffusion timescale. However, the timescale for 
the disk to get close to the steady state is much longer than the 
viscous timescale. In order to explicitly show this, we plot, 
in Figure [2] fQs(r), defined as the time the surface density at 
radius r takes to reach 95% of its limit value (in solid line). 
Here, QS stands for quasi-steady state. We compare those 
timescales to the viscous timescale, 

tew = r 2 /v , (3) 

which is a timescale for how quickly diffusion can change the 
density at a particular radius. Here, v is the kinematic vis- 
cosity of the gas, for which we adopt th e turbulent (or eddy) 
viscosity of IShakura & Sunvaevl (11973b . v = ac s H, with c s 
being the sound speed and H the disk scale height at a given 
radius. The viscous diffusion timescales are plotted in Fig- 
ure [2] as the dashed lines. The viscous diffusion timescales 
are much shorter than ?qs and are not, therefore, represen- 
tative of the time the disk takes to settle in a quasi-steady 
steady after a long period of constant injection rate. For in- 
stance, a disk with a = 1.0 takes 30 years to reach 95% of 
the limit density at 107?*, whereas the diffusion time for this 
radius is only 1.7 years. Although most of the changes oc- 
cur at a viscous diffusion timescale, it follows from the above 
that the disk requires a longer time to be filled up and stabi- 
lize in its outermost reaches. It is probably difficult to find 
a Be system in which the entire disk is close to steady state, 
since this would require an exceptionally long and stable pe- 
riod of mass loss from the central star. However, ( Tauri and 1 
Del could probably be emblematic cases of this behaviour and 
more generally, late-type Be stars, since they are supposedly 
more st able than early-type ones as seen from photometric 
studies dCuypers et al| [l989t iBalona et al.lfl992t lsiagg|[T987t 
Uones et al.ll20TlF We conclude that the timescales for the 
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disk to fully redistribute the material ejected from the star, 
(§0, are much larger than the viscous d iffusion timescales. 

From observat ions (e.g. .[Wat ers 1986) and theoretical con- 
siderations (e.g. jBiorkman & Carci ofi 2005), Be disk models 
often use a power-law to describe the radial mass density pro- 
file of the disk, i.e., 

p(r) oc r~" . 

Let us now discuss what are the model predictions for the 
value of n. Assuming vertical hydrostatic equilibrium and an 
isothermal gas, the mass density is related to the surface den- 
sity by 



P(r, z) 



Z(r) 



y/2nH(r) 



exp 



2H 2 (r) 



(4) 



where r and z are the usual cylindrical coordinates. Since 
the isothermal scale height increases as r 1 5 and the steady- 
state surface density decreases as r~ 2 , we obtain the famil- 
iar result that the ste a dy sta te volume density falls as r -3-5 . 
ICarciofi & Biorkmanl ( 120081) investigated how the disk tem- 



perature structure affects the viscous diffusion. The combi- 
nation of the radial temperature structure, disc scale height, 
and viscous transport produces a complex radial dependence 
for the disk density that departs very much from the simple 
n — 3.5 power-law. This is certainly a source of explanation 
to the wide range of n values reported in the literature. But 
in the present paper, our results suggest that there is an addi- 
tional, alternative explanation for this. 

In Figure |3] we show the local index of the surface density 
profile defined as m = — dla(Z)/dln(r). From Eq. (|4]i we see 
that local index of the volume density, «, is related to m by 
n — hi + 1.5, for isothermal disks. A large value of m at a par- 
ticular location indicates that the disk is far from steady state 
at that location. As the disk builds up, the slope progressively 
reaches the limit value of 2, beginning at the inner disk and 
spreading out to larger radii. The time required for the outer 
disk to reach the m = 2 slope can be quite long (decades), 
even for large values of a. Therefore, if one estimates the 
slope of the density profile (for instance, via the IR SED), one 
is likely to find values of m > 2 (or, equivalently, n > 3.5). 
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Fig. 2. — Time required for the disk density to reach 95% of the quasi- 
steady state value (fQs (solid lines) and viscous timescale, t^g (dashed lines), 
as functions of the distance from the star and computed for different values 
of or: 0.1,0.3, 0.5, 0.7 and 1.0. 

This could therefore account for the wide range of n values 
in the literature. Note that this is only true for the disk build- 
up phase since, as presented below, the density slope for the 
dissipation phase is quite different. 

3.1.2. Disk dissipation 

When mass injection from the star ceases, the angular mo- 
mentum supply to the disk stops and the inner disk starts to re- 
accrete back onto the star. Reaccretion of material from the 
inner disk occurs because the turbulent viscosity transports 
angular momentum outward, transferring angular momentum 
from the inner to the outer parts of the disk. 

To study the disk dissipation we start from a pre-existing 
disk and turn off mass input from the star. The temporal evo- 
lution of the surface density is shown in the right panels of 
Figure Q] for two values of a (0.1 and 1). As soon as the 
mass input is turned off, the inner disk re-accretes back onto 
the star, giving rise to a stagnation point , i.e. a region in the 
disk where the radial velocity is zero. The stagnation point 
slowly propagates outwards and the density of the inner part 
decreases in a homologous way with time, i.e. the density 
structure within the stagnation point remains self-similar dur- 
ing this phase. The speed in which the stagnation point propa- 
gates outwards is shown in Figure|4l it is, as expected, strictly 
proportional to a. For instance, for a — 1 the stagnation point 
reaches about 100 /?* in 1 year. This time, at which accretion 
starts at 100 is actually shorter than the viscous time (~ 
5.5 years). 

Usin g the formalism presented in iBiorkman & Carciofil 
(120051) . it is possible to demonstrate that the surface density 
of an accretion disk in the steady-state regime is given by 



E(r) = 



MV cnt fi, 




(5) 



where c s is the sound speed. 

In Figure [5] the shape of the simulated density profiles for 
a decaying disk is perfectly matched by the analytical curve 
of the surface density for a steady accretion state using /?o = 



0.85/?*. This parameter is particularly important because it 
determines the shape of the density in the inner part of the 
disk. As time goes on, the stagnation point goes outwards, 
the accretion region grows and a larger fraction of the disk 
assumes a density distribution similar to Eq. For example, 
15 years after accretion starts, the first 50 are already in a 
steady accretion state. 

We conclude that during the disk dissipation the disk as- 
sumes a dual behavior, with two distinct regions (inner accre- 
tion and outer decretion) separated by the stagnation point. In 
the decreting part, the index m of the surface density is about 
2, whereas in the accreting part the slope goes from 1.5, which 
is the value for a steady-state accreting disk, to negative val- 
ues closer to the star (see Figure[6j. 

3.2. Periodic injection rate 

Having studied the limiting case where the T; n » tj, we 
now proceed to study cases for which Tj n ~ t^. To study how a 
disk forms and decays under the action of a variable mass loss, 
we studied some test-cases of periodic injection rates. In those 
models, we successively turned the injection rate on and off 
for different periods of time and computed the corresponding 
surface densities. Some examples are shown in Figure [7j For 
such Ti n ~ Td cases, we expect that the injection rate variation 
will generate a combination of the two previously described 
fundamental cases of disk build-up and decay. To investigate 
those periodic scenarios we used three parameters: the period, 
the a viscosity parameter and the duty cycle (hereafter DC, 
the time spent in an active state as a fraction of the total time). 

On Figure [7j we can see the effect on the surface density of 
a periodic material supply in a disk for different DCs. This 
succession of decay and build-up phases results in a com- 
plex temporal behaviour of the density profiles. For the same 
epochs, the curves are very specific for each DC because they 
depend on the previous dynamical history of the disk. How- 
ever and as a general statement, we can say that for short 
timescales (i.e., within a given cycle) there is again a two- 
component structure: 



• the inner part of the disk that oscillates ; its surface den- 
sity is consequently flapping up and down every period 
in fast reaction to the injection rate variations, and 

• an outer part of the disk that is less affected by this in- 
termittent mass injection but undergoes delayed, much 
slower variations. 

In order to illustrate the long-term changes of the disk, it is 
interesting to look at the surface densities computed for the 
same phase along different cycles. Figure [8] shows the sur- 
face density for four different phases (0.05, 0.45, 0.55 and 
0.95) at several cycles. In general, cycle to cycle changes of 
the inner disk (say, within 10 /?*) are small. This is to be 
expected because the inner part responds quickly to changes 
of the mass injection rates. The outer part, however shows a 
steady growth in density as the disk is slowly filled up by the 
mass supplied in the previous cycles. 

Given that the density is slowly but steadily increasing in 
the outer part, one can wonder whether the density of a disk 
undergoing a periodic injection rate will ever reach a limit 
value in its outermost reaches. Figure [9] shows the temporal 
variation of the surface density at three different radii in order 
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Fig. 3. — Index m of the surface density, Ear "', as a function of the radius during disk build-up phase (left: a = 0.1, right: a = 1.0). Epochs from top to 
bottom are: 0.1, 0.5, 1.0 , 2.5, 5.0, 15.0 and 50 years. See Figure 1 and text. 
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Fig. 4. — Position of the stagnation point, i.e. size of the accretion region, 
as a function of time for a dissipating disk for different a (bottom to top: 0. 1 , 
0.3, 0.5, 0.7 and 1.0). 



to answer to this question. It is demonstrated that indeed the 
disk reaches a limiting value which is DC x l,o(R+/r) 2 . 

It is interesting to note the density oscillation at 3 R* in Fig- 
ure [9] a behavior that is not observed at larger radii. This 
happens because for small radii the disk alternates between 
decretion and accretion, whereas the outer disk never expe- 
riences accretion, only decretion. This is quantified in Fig- 
ure [10] which shows the position of the stagnation point as a 
function of time for three values of a. 

As mass injection stops, the accretion region increases in 
size until mass injection starts again. The higher the a values, 
the higher the viscous torque, the further the accretion region 
expands. The maximum extension of the accretion region sta- 
bilizes over time to reach 40 R+ (a = 1 .0) after ~ 40 years. To 
summarize, we find that long periods, low duty cycles (equiv- 
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Fig. 5. — Surface density profiles from singlebe simulations (solid lines) 
and from the analytical expression of Equation \5\ (dashed line) in the disk 
dissipation scenario. The thick line represents the quasi-steady state density 
profile of Eq. (2) and the black dot indicates the stagnation point, a = 0.5 
and the four epochs shown are from top to bottom: 1, 5, 15 and 35 years. 

alent to short time outbursts) and high a values are needed to 
build the largest accretion regions. 

In this subsection we saw how a disk builds in reaction to a 
periodic scenario. The evolution of the surface density radial 
profile depends strongly on the periodicity parameters (DC 
and period) and also on a. Contrary to the simpler scenarios 
explored in Sect. 13.11 here the surface density depends on the 
details of the previous dynamical history of the system. 

3.3. Episodic mass injections 

Even if periodicity is observed in lightcurves, they are also 
sometimes just random dMennicken t et al. 2002). As an ad- 
ditional class of dynamical scenarios for which r m ~ Td, we 
investigate cases showing of a sudden increase in the mass 
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Fig. 6. — Same as Figure fJ]for a dissipating disk. Shown are results for two values of a (left: or= 0.1, right: a =1.0) and for the following epochs (top to 
bottom): 0, 0.1, 0.5, 1., 2.5, 5, 15 and 50 years. The horizontal and vertical dashed lines mark the final position of the maximum surface density. 



injection rate (i.e. an outburst). 

A representative model is shown in FigureQj] Here, a con- 
stant mass injection for 20 years is followed by an outburst 
with twice the mass injection rate lasting 0.2 year. After the 
outburst, the injection rate was set to zero. The surface density 
during the outburst reveals an important property of outburst- 
like scenarios: the outburst affects mostly the inner part of the 
disk, and how far out the surface density is affected depends 
on the length and strength of the outburst, as well as on the a 
parameter. 



4. PREDICTIONS OF THE DYNAMICAL PROPERTIES OF THE 
PHOTOMETRIC OBSERVABLES 

In the previous section, we described the evolution of the 
disk structure for different dynamical scenarios. This part of 
this paper is devoted to hdust predictions of the variability of 
photometric observables. A large variety of observables can 
be derived from the hdust simulations. To structure the analy- 
sis, we choose to present in this paper only the most important 
photometric results. Polarimetric, spectroscopic and interfer- 
ometric results will be analyzed in future publications. 

Photometry is a powerful tool to use when observing vari- 
able stars because fluxes at different wavelength s allow prob- 
ing th e disk at different radii (e.g., Figure 6 of Ide Wit et alJ 
120061) . As we want to follow the outward growth and dis- 
sipation of the disk, it is critical to understand from which 
part of the disk comes most of the emitted fluxes at different 
bands. Figure [T21 shows the normalized radial flux distribu- 
tion of a disk around a Be star at various wavelengths. At a 
given band, the intersection of the flux curve with the horizon- 
tal dashed line marks the position in the disk whence about 
95% of the continuum excess comes. For instance, we see 
that the V-band excess is formed very close to the star, within 
about 2 stellar radii, whereas the excess at 1 mm originates 
from a much larger area of the disk. To understand this plot, 
the disk can be regarded as a pseudo-photosphere around the 
star, whose effective size increases with wavelength. Note that 
these predictions for the formation loci are for an inclination 
angle i = 30°. As we will see below, i bears important effects 
on the emergent flux. Finally, Figure[T2]is in broad agreement 
with interferometric measure ments of disk sizes for near-IR 
and mid-IR wavelengths (e.g. lGies et alj20O7tlMeilland et alJ 



120091) . However, the extent of the A^-band continuum emit- 
ting regi on seems system a tically smaller than expected for ex- 
ample in lChesneau et alJ (120051) (see Figures 6 and 7), but in 
agreement with observations. 

Based on the previously described dynamical scenarios, we 
computed various photometric observables at different wave- 
lengths with hdust in order to study the mass redistribution 
process at different locations in the disk. In the following, we 
selected some results for the V-band, /iT-band and millimetric 
domains. 

4.1. V-band photometry 

The V-band magnitude in the disk is controlled by two pro- 
cesses: gas emission and absorption of the stellar radiation. 
From Figure [12] we expect that only the densest inner part of 
the disk (within r < 2R+) will affect the flux at this wave- 
length. Also, these two mechanisms will have a different im- 
pact on the observables depending on the inclination angle. 
At low inclination angles, the disk is seen face-on, and lit- 
tle absorption of the stellar radiation is expected. At high 
inclination angles, the disk is seen edge-on, and absorption 
plays a stronger role. Figure [T3l shows examples of V-band 
light curves com puted f or the dynamical scenarios described 
in § 13.1.21 and § 13. 1 . 1 1 for three different inclination angles. 
We plot AV, the difference between the total flux and the pho- 
tospheric flux in the V band. 

For the build-up case, we can check the above statements on 
the inclination angle. For i = 0°, for instance, the disk devel- 
ops a strong V excess of -0.4 magnitude. For i = 90°, the disk 
causes a flux decrement of about 0.2 magnitude, as a result 
of absorption of photospheric light. A balance between emis- 
sion and absorption is met for an inclination angle of about 
70°. This fact is wavelength-dependent as we will see that the 
balance is met at different values of i for other spectral bands. 

For a given i, the timesc ale f or the V evolution depends 
strongly on a. As seen in §13.11 the higher the a, the faster 
the disk will reach the limit density, so the faster the V-band 
variations will be. We should state here that the asymptotic 
values reached by AV depend on the density in the inner part 
of the disk and thus on the particular choice of So- Table [3] 
summarizes how long AV takes to reach 95% of its asymptotic 
value depending on a and the inclination angle of the disk. 
As V is mostly affected by the inner disk, those values are not 
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Fig. 7. — Left: Periodic mass injection rate with a 2 year period for different DCs (top: 10%, middle: 50%, bottom: 90%). Right: Corresponding surface density 
profile for o-=0.5 at chosen epochs. Bottom to top: 5.9, 6.01, 6.05 and 6.1 years. Black dots indicates the stagnation point. 
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Fig. 8. — Surface density profile for a periodic mass injection rate with a 2 year period, a=0.5 and DC = 50%. The four panels show four different phases at 
6 different cycles: 1st, 2nd, 3rd, 6th, 26th and 51st. Top left: phase 0.05, top right: phase 0.45, bottom left: phase 0.55, bottom right: phase 0.95. Black dots 
indicates the stagnation point. 



sensitive to different disk sizes as long as Rj > 2/?*. 

Figure Q~3] illustrates another important prediction of the 
VDDM: the timescales for disk growth are shorter than the 
timescales for disk dissipation. This is due to the fact that 
while at disk growth, the timescales are set by the matter re- 
distribution within a couple of stellar radii only. At disk dissi- 
pation, the timescales are controlled by the reaccretion from a 
much larger area of the disk. This prediction is confirmed by 
observations dCarciofi et al. 2 0121) . 

Figure [14] shows the predictions for the scenarios with pe- 
riodic mass injections. Contrary to the case of a steady disk 
growth in which the inner disk forms without perturbation and 
steadily reaches the limiting density, the inner disk has not the 
time to fully develop in the case of periodic scenarios. This 
results in lower amplitudes for the variations in AV. Clearly, 
this effect is more pronounced from smaller DCs. 

The diverse morphology of the V-band lightcurves shows 
that they are quite specific of the mass decretion scenario, 
in this case the DC and cycle length. They are also quite 



TABLE 3 

Time in years required for V to reach 95% of its asymptotic 

VALUE. 



a I Angle (deg) 30 70 80 85 90 



0.1 


27.8 


29.3 


<0.1 


26.0 


71.0 


90.0 


0.3 


11.4 


12.7 


<0.1 


9.7 


23.6 


24.9 


0.5 


7.1 


7.8 


< 0.1 


7.9 


12.0 


14.0 


0.7 


4.1 


4.3 


< 0.1 


4.4 


9.7 


11.4 


1.0 


2.7 


2.8 


<0.1 


3.4 


7.1 


8.1 



dependent on the a parameter. Therefore, follow-up obser- 
vations of the V-band magnitude represent a powerful tool 
to infer these paramete rs. Its diagno s tic po tential was re- 
cently demonstrated by ICarciofi et alj d2012l) who used the 
lightcurve of the Be star 28 CMa to measure, for the first time, 
the viscosity parameter in a Be star. It is important to keep in 
mind, however, that V-band variations are rather insensitive 
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Fig. 9. — Temporal evolution of the surface density at three different disk 
locations: 3 (top in solid line), 10 i?* (middle in dotted line) and 100 R+ 
(bottom in dash-dotted line). The corresponding periodic mass injection rate 
has a 2 year period, o>=0.5 and DC = 50%. The horizontal lines correspond 
to the limit value presented in $ 13.1.11 The horizontal lines correspond to the 
limit value at that location (Eo * i" ) times the DC. 




Fig. 10. — Temporal evolution of the accretion region size for different val- 
ues of a (top to bottom: 1.0, 0.5 and 0.1). Shown are results for a period of 2 
years and a DC of 50%. 



to the outer disk parts. 

If lightcurves are of great importance to study Be stars, 
color-color and color-magnitude diagrams can also bring fun- 
damental elements to the analysis because they are character- 
istic of different radiative transfer effects occurring in the disk 
(absorption, emission, etc) at different radii. The disktem per- 
ature is typically 60% of T e g (ICarciofi & Bior kman 20061), so 
the spectral shape of the disk emission is redder than the pho- 
tospheric spectrum. Therefore, colors are a good tracer of the 
amount of stellar radiation that is reprocessed towards longer 




Fig. 11. — Top: Mass injection history of an outburst-like model. Bottom: 
Corresponding surface density profile for representative epochs, as indicated. 
For this model, a = 0.5. 



TABLE 4 

Time in years required for K to reach 95% of its asymptotic value.. 



a 1 Angle (deg) 





30 


70 


80 


85 


90 


0.1 


26.3 


28.5 


92. 


> 100. 


> 100. 


<0.1 


0.3 


9.8 


11.2 


24.2 


33.5 


95.0 


<0.1 


0.5 


6.0 


6.7 


16.6 


31.3 


46.6 


<0.1 


0.7 


4.0 


4.4 


9.9 


20.2 


27.6 


<0.1 


1.0 


2.9 


3.3 


8.4 


12.1 


19.7 


<0.1 



wavelengths by the disk. 

In Figure we show examples of U-B/B-V color-color 
and B-V/V color-magnitude diagrams, illustrative of disk 
growth and dissipation. Starting from zone A (no-disk and 
a small color decrement) the system follows a reddening path 
towards an asymptotic value (located in zone B). When mass 
loss is turned off, the system follows a blueing path back to 
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Fig. 12. — The formation loci of different observables in the continuum 
emission. The calculations assume the parameters listed in Tables 1 and 2. 
The results were derived for an inclination angle of 30 °. Plotted is the ratio 
between the observed flux to the maximum flux F max , defined as the flux of a 
model with a disk outer radius of 1000 R+, as a function of the distance from 
the star. 
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Fig. 13. — V-band lightcurves associated with disk build-up (top panel) and 
dissipation (bottom panel). The dotted, dashed and dot-dashed lines represent 
the lightcurves for inclination angles of 0° (face-on), 70° and 90° (edge-on), 
respectively. The black, red and blue colors represent models for ci'=0. 1, 0.5, 
and 1.0, respectively. The solid black lines indicate the asymptotic value of 
AV. 



A that is different than the previous path , therefore forming a 
loop in the color-color diagram. The different paths from A to 
B and from B to A come from the fact, discussed above, that 
the disk grows at a faster rate than it dissipates. 
The building up and decay phases (and by extension pe- 
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Fig. 14. — Same as Figure [T5] for the periodic scenarios of Figure [7] (top: 
DC=10%, middle: DC=50%, bottom: DC=90%). 



TABLE 5 

Time in years required for the mm magnitude to reach 95% of 
ITS asymptotic value.. 



a 1 Angle (deg) 





30 


70 


80 


85 


90 


0.1 


26.2 


26.7 


28.4 


33.4 


38.3 


40.1 


0.3 


9.4 


9.3 


9.7 


12.9 


15.4 


16.6 


0.5 


5.8 


4.6 


6.4 


7.4 


8.1 


8.4 


0.7 


3.9 


2.6 


4.2 


5.3 


6.2 


6.5 


1.0 


2.5 


2.7 


2.9 


3.5 


3.8 


4.0 



riodic injection rates) show clear loop-like signatures in the 
color-color U-B/B-V and color-magnitude diagrams. As 
the morphology of those loops depends strongly on the incli- 
nation angle, those diagrams represent a tool to estimate this 
parameter on top of the decretion history at play. Photon noise 
is particularly visible on the 90 degree figures since the pho- 
tometric signal is much lower than for the other inclination 
angles where absorption plays a lesser role. 

4.2. K-band photometry 

The excess flux in the IR comes mostly from a part located 
up to several stellar radii from the star (~ 6R* in /if -band, 
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Fig. 15. — Top panels : U-B/B-V diagrams for a 50 year long disk build-up followed by a 50 year long dissipation. Each panel shows the results for a 
different inclination angle. Epochs are indicated in years. The path A to B represents the build up phase and the B to A the dissipation phase. Bottom panels : 
corresponding AV/B-V diagrams. 



see Figure [T2l. This is also an interesting region because this 
is where we observe the h ighe st density variations for a peri- 
odic decretion scenario (§ 13.21 ). so we expect clear signatures 
in the /T-band lightcurves. We mention here that a variability 
of a different sort in the disk structure is to be expected due 
to thermal effects that are n ot simulated in this work. T hose 
effects have been studied in lCarciofi & Biorkmanl ( 120081) and 
Uones etalJd2008l) . Chan ges in d ensity can be important in th e 
inner disk (e.g. see Figure 4 in ICarciofi & Biorkmanll2008h . 
Their consequences upon short-wavelength observables re- 
main to be investigated and constitute an important perspec- 
tive for this work. 

Figure [16] shows the influence of the inclination angle and 
the a parameter on the /T-band flux for the build-up and dis- 
sipation phases. Examples of lightcurves for the periodic sce- 
narios are shown in Figure [17] While the behavior is similar 
to the V-band case, there are three important differences: 

1 . The asymptotic value of the excess reached for i = 0° 
is much larger than for the V-band. This is simply a 
result of the fact that the size of the pseudo-photosphere 
grows with wavelength. 

2. Both the growth and decline rates of the /if -band 
lightcurves are slower than in the V-band. 

3. The inclination for which emission and absorption bal- 
ance each other is close to 90°. This means that, for sys- 
tems seen at that such an inclination, even if the disk is 
experiencing long growth and dissipation phases, they 
will be barely detectable in this band. 



From the above, we see that the combination of V and K 
bands has a strong diagnostic potential. For instance, it prob- 
ably allows for constraints to be put on the inclination an- 
gle. More importantly, the fact that each lightcurve is a diag- 
nostic of the mass redistribution timescale at different radii, 
combining V and K lightcurves opens the intriguing possibil- 
ity of detecting possible variations of the a parameter with 
distance from the star. For example, a mismatch between 
the growth and/or decline rates of the model and observed 
lightcurve could be indicative of a radial variation of a. 

In the dissipation lightcurves, an interesting phenomenon 
happens for i = 90° (lower panel of Figure [Toll: as the inner 
disk empties, the /iT-band emission of the disk decreases and, 
at the same time, the outer parts still absorb radiation both 
from the star and from the inner disk. Consequently, AK be- 
comes positive and then slowly goes to zero as the outer disk 
disappears. 

4.3. Millimetric photometry 

The millimetric (mm) emission (at A = 1mm) comes from 
an extended outer part of the disk (mainly between 5 and 
50R*). This implies that the mm lightcurve is highly sensi- 
tive to the disk size. We thus explored four disk sizes in the 
range 5 — 100/?*, to study the effect of disk truncation by 
a secondary star. Figure [18] shows the mm lightcurves during 
the build-up phase of these disks of different sizes. On this 
figure, we can see that the time required for the mm magni- 
tude to reach a stable value depends strongly on the disk size. 
The bigger the disk, the higher that value. 

The magnitude variations at 1 mm computed for the build- 
up and decay phases are plotted in Figure [191 

As in the V and K bands, the mm emission is produced 
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Fig. 16. — Upper two panels: same as Figure [T3l but for the K band. Lower 
panel: a zoom into the dissipation curve allows to see that AA^ becomes posi- 
tive for i=90° . 
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Fig. 17. — ff-ba nd lightcurve for the periodic scenario. Legend is the same 
as for Figure [T3l 



quite fast after the mass decretion started. Table [5] shows how 
fast the mm disk excess reaches 95% of its limit value. As 
a major difference with the V and K band, during the decay 
phase, it takes a very long time before the emission decreases 
significantly: in 50 years, the excess at "decreases by 1.5 
magnitude for a = 0. 1 . 

Figure|20]shows the mm lightcurves during a periodic mass 
injection scenario for different angles and a parameters. The 
combined influence of the inclination angle, DC and a on 
the mm lightcurves is not trivial. However, Figure [20] clearly 
shows that for low a values, the periodic structure imposed by 
the decretion scenario is not visible in the lightcurves for high 
DCs (middle and lower panels of Figure [20l. The reason is 
that the area responsible for the mm emission is so large that 
if the disk starts to vanish slowly from its inner parts for a few 
months and then gets filled up again, the resulting lack of den- 
sity will propagate but will remain unnoticeable for the over- 
all mm emitting area. However for low DCs (upper panel of 
Figure [20l, the effect of a periodic injection rate is visible be- 
cause the mass supply injected into the disk is comparable to 
the overall mass of the mm emitting region of the disk. How- 
ever, this effect is more obvious for high a values where the 
density variations are more important because of large vari- 
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Fig. 18. — Lightcurves at the mm domain for a build-up phase under a con- 
stant injection rate (i = 80°, a = 0.1). The black, red, blue and green colors 
respectively represent the lightcurves for different disk sizes: 5, 10, 20 and 
100 Sy, respectively. 



ations of the accretion region size (see discussion in § 13.21 ). 
Another interesting fact that we can observe on this figure is 
the temporal delay between the injection rate variation (turns 
on at 8 years) and the change in the lightcurve shape. Depend- 
ing on the a value, the change of the disk state is visible in the 
lightcurves after 1 (a =1.0) or 2 (a = 0.5) months. 
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Fig. 19. — Same as Figure fPH but for a wavelength of 1 mm. 

We can summarize then that the mm lightcurve is a good 
observing tool to infer a disk size but not to determine a short- 
term mass injection history neither to estimate unambiguously 
the inclination angle and the a parameter. 

4.4. Photometric consequences for episodic models 

The lightcurves corresponding to the episodic model are 
presented in Figure [2T1 

The episodic model involves a 2 month outburst of stellar 
mass loss whose mass injection rate is twice the one the disk 
underwent during its 20 first years of building. The resulting 
lightcurves of such an event can be explained with the same 
arguments as previously. It is interesting to compare the ab- 
solute excess values due to the outburst with the excess at the 
end of the previous build-up phase. The biggest difference is 
in the V-band where we see that an increase of the injection 
rate by a factor 2 results in a maximum magnitude value of 
-0.55 at 0°. The outburst is a little less conspicuous in the K- 
band, and barely visible in the mm lightcurves. This illustrates 
the important fact that the best indicator of the real time mass 
injection rate is photometry at short wavelengths. 

5. DISCUSSION 
5.1. Disk structure 



As outlined in § 13.1.11 the density structure of the disk is 
often parametrized in the form of a power law (p(r) oc r _n ). 
For a well-developed disk, with a flaring parameter of 1.5, the 
expected value is n = 3.5. From visible interferometric ob- 
servations of x Oph, iTvcner et alj d2008l) foun d a best-fitting 
value of n = 2.5, whereas Ijones et al.l (120081) reports values 
of 4.2, 2.1 an d 4.0 for k Dra, B Psc and v Cyg, respectively. 
IWatersi (119861) report values in the range 2 - 3.5, based on an 
analysis of the IR SED of several Be stars. Thus a wide range 
of values for n is seen in the literature. 

We found (ICarciofi & Biorkmanll2008l) that non-isothermal 
viscous effects can be one possible explanation for large scat- 
ter of the density radial exponent reported in the literature. 
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20. — Mm- band lightcurves for the periodic scenario. Legend is the 
as for Figure FBI 



However, another, and, in view of the perpetual variability of 
Be star disks, quite an attractive mechanism lies in the steep- 
ening of the index in decretion phases (« > 3.5 in case de- 
cretion starts with an empty disk) and flattening of the index 
in re-accretion phases (n < 3.5 in case decretion starts with a 
fully developed disk). Moreover, the index changes with dis- 
tance from the star, and thus any determination of the index 
will be sensitive not only to the decretion history, but as well 
to the wavelength for which it is determined. 

5.2. Observed Galactic Be stars 

Data for some example stars were taken from published 
photometric databases. The three best examples found were 
28 CMa, representing a low i Be star with a decaying disk, 
88 Her, a shell star (i.e. at high ;) which has shown signs of 
two about one-year long mass transfer events during the al- 
most ten years of observations, and s Cap, another shell star, 
which alternated between a bright and faint state every other 
year for about eight years of observation. For the discussion, 
we note that colors are bluer when they are more negative. 
Data are shown in Figure [22] 

The dataset of the almost pole-on Be star 28 CMa has al- 
ready been described by IStefl et ail ((2003). The observa- 
tions started in a phase of mass transfer, reaching a state of 
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Fig. 21. — V-band, K-bwd an d m m lightcurves for the episodic scenario. 
Legend is the same as for Figure [T3l 



a well developed disk around MJD=46000 in which B-V is 
reddest and U-B is bluest. Then the disk started decaying 
about a year later, i.e. mass transfer weakened or ceased. 
The color-color path during decreasing brightness is bluer 
than during increasing brightness, to finally reach a state in 
which B-V is bl uest and U-B is reddest. This agr ees with 
lHarmane"c1 (1 1 9831) . as well as with the results of § 14. H and Fig- 
ure [15] Another star showing such behaviour of secular color 
changes during d isk built-up and decay is KDra in the 1980s 
dJuzaet al.lll994h . 

The shell stars 88 Her and eCap show a differ e nt be- 
havi our. The data were tak en from iPavlovski et alj (119971) 
and iMennickent et alJ (119941) . respectively, and if necessary 
converted to the UBV sy stem using the relations given by 
lHarmanec & Bo zic (2001). These are bluer in U-B when 
having little circumstellar material (meaning when they are 
bright due to lack of absorbing gas). As material is decreted, 
J/-Bclearly and B-V tendentially become redder. The U-B 
behaviour is again in good agreement with § 14.11 and Fig- 
ure Q3J while B-V is not. However, the amplitudes, both ob- 
served and predicted, in B — V are much smaller than in U-B. 
As mass injection ceases, the star 88 Her returns to the bright 
state on a bluer path, for s Cap this cannot be said with cer- 
tainty. Also this is in agreement with the modeling. 



Most other stars in the photometric databases show a gen- 
erally less clean picture, with the notable exception of 48 Lib. 
In that star U-B and B-V are positively correlated, which 
is not seen in any of the model co mputations. We no te, that 
such a correlation is also given by lHarmane"c1 (119831) as the 
normal behaviour of shell stars. However, the colors not only 
correlate well with each other, but as well with the state of 
the long-term V/R va riation (see e.g. IMennickent et alJI 19981 : 
iMcDavid et al.ll2000l) . while at the same time the Ha equiv- 
alent width and emiss ion height was constant (Figure 1 1 of 
lHanuschik et al.|[1995[) . This suggests that the color and mag- 
nitude changes observed in 48 Lib are not linked to variability 
of the mass injection rate (and hence the radial density struc- 
ture), but instead to the azimuthal structure of the disk gov- 
erned by the one-armed density wave, seen at different aspect 
angles over the years. Since the B-V amplitude of 48 Lib is 
twice as high as the ones of 88 Her or eCap, it is possible 
that a correlation as given by Harmanec is governed rather by 
V/R-cycle related changes than by changes related to the mass 
injection. 

5.3. Comparison with previous models 

Already in 119781 iPoeckert & Marlboroughl published 
theroretical computations for brightness and color variations 
for a set of steady-state disk models, i.e. assuming the stable 
equilibrium of a fully developed disk. Even if the basic 
assumption of those models have been revised considerably 
in the more than 30 years since, their findings concerning the 
photometric behaviour remained the only available parameter 
study until now (their Figures 29 to 34). In quantitative terms, 
their findings do not differ strongly from the ones presented 
h ere for disks after a long built-up phase, i.e. as described in 

IPoeckert & Marlboroughl also find an inclination at which 
the brightness change is zero. The inclination depends on 
wavelength, but the value at which the cancellation of ab- 
sorption vs. addi tional emission occurs is generally smaller. 
For the V-band, IPoeckert & Marlboroughl find an invariant 
magnitude for an inclination of about 45 to 60 °, while 
here it is 70. At a wavelength of 2.7pm, an exactly 
equator-on disk still redu ces the brightness in the model of 
IPoeckert & Marlboroughl while here it is found that long- 
wards of about the /T-band even a perfectly equator-on seen 
disk will increase the brightness of the star. Only at 10pm 
both model predictions agree again qualitatively for an edge 
on disk, in that they both brighten. Another importa nt quanti- 
tative difference is that for lPoeckert & Marlboroughl the (pos- 
itive) magnitude difference in the edge-on case is about twice 
the value of the (negative) magnitude difference when seen 
pole-on (AV et i ge - on = -2AV po i e - on ), whereas this ratio turns 
out to be reversed in this work (AV ei jge-on = -0.5AV /TO / f _ OJ1 ). 

In ter ms of color-color variatio ns, for a comparable base 
density, Poe ckert & Marlboroughl derive somewhat smaller 
changes than here, but it is important to note that these 
changes have the same sense: 

• For a disk seen face-on, V is brighter, U-B is bluer, 
while B-V becomes redder (vs. a diskless star). 

• For a disk seen edge-on, V is fainter, U-B is redder, 
while B-V becomes bluer (vs. a diskless star). 

As already pointed out above, the second point is at vari- 
ance with the observational correlation given by lHarmanecl 
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Fig. 22. — Color and magnitude variations for three Be stars. The filled circles are bimonthly (28 CMa) or annual (88 Her, sCap) means. The line sketches in 
the respective uppermost panels indicate the direction of the path in time. 



dl983l) . which, for stars with "inverse correlation between H I 
emission and luminosity" states that the "fading in V is ac- 
companied by the redden ing of both U—B and B-V". Both 
iPoeckert & Marlborough! and this work find a blueing of B— V 
if the fading in V is caused by an axisymmetric increase of 
disk density and mass. This lends some support to the hy- 
pothesis that the observed correlation might in fact not (fully) 
be governed by growth or decay of the disk. 

The quantita tive differences are most lik ely because the disk 
considered by IPoeckert & Marlborough! was quite a simpli- 
fied one compared to the current understanding of circumstel- 
lar disks, e.g. it included no flaring and was isothermal (and 
hotter) than the disk considered here. In qualitative terms, 
however, we find good agreement with the previous results. 

6. CONCLUSIONS 

The goal of this paper is to report on photometric predic- 
tions derived from the modeling of viscous decretion disks 
that undergo variable mass injection rates. We thus simulated 
the disk evolution under the influence of different dynamical 
scenarios. In order to understand the evolution of the disk 
structure in details, we looked at the temporal variations of 
some fundamental quantities derived from the surface density 
profiles. There are several timescales at play in the disk that 
helps to understand its evolution. To provide the reader with 
some reference numerical values, we made a quantified com- 
parison between the viscous diffusion time and the time the 
disk surface density requires to reach its limit value. Coming 
also from the study of density profiles, we brought some the- 
oretical arguments to explain the variety of power-law index 
values reported in the literature. We then presented V-band, 
/if -band and mm lightcurves at some epochs of the dynam- 



ical scenarios. The constant, periodic or episodic decretion 
scenarios generated characteristic variations on those photo- 
metric observables. We noted that the general behaviour of 
the V and K excesses is quite similar. These lightcurves vary 
accordingly to the inclination angle in a first extent, then to 
the a parameter. They consequently represent a good tool to 
estimate those parameters. Moreover, as they result from an 
emission produced in the first stellar radii, they are therefore 
a good way to infer, at a several days timescale, the mass loss 
history of the central star. We found that mm lightcurves are 
more useful for disk size determination. However, there is 
some degeneracy, i.e. several sets of decretion scenarios and 
parameters could describe a same given short-term observed 
lightcurve. Consequently we stress that the more observables 
with a long time coverage, the easier to infer a dynamical sce- 
nario and physical parameters of the system. We finally com- 
pared our results to reported observations and results from 
previous modelings. In order to further test the o--disk the- 
ory, we plan to compare our models with lightcurves from a 
representative sample of Be stars observed at different wave- 
lengths. Moreover, additional predictions on spectroscopic, 
polarimetric and interferometric observables will be presented 
in further publications. 
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